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Abstract 

We present next-to-leading order QCD corrections to production of two W bosons at the 
LHC in the Randall- Sundrum model. Various kinematical distributions are obtained to 
order in QCD by taking into account all the parton level subprocesses. We estimate the 
impact of the QCD corrections on various observables and find that they are significant. 
We also show the reduction in factorization scale uncertainty when 0{as) effects are 
included. 
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The main aim of the upcoming Large Hadron Colhder (LHC) is the search of the 
missing piece of the standard model (SM) ie. the Higgs boson and the existence of 
new physics which offers the solution to the hierarchy problem of SM. In this direction, 
there exist many models based on ideas of supersymmetry (SUSY), extra dimensions, 
technicolor etc. The possible existence of new spatial dimensions beyond 3+1, came 
from early works of Kaluza and Klein in which they postulated a fifth dimension to 
unify electromagnetism and gravity [1], but the renaissance of extra dimensions began 
with the proposals of Arkani-Hamed, Dimopoulos and Dvali (ADD) [2], and Randall and 
Sundrum (RS) ||3]. 

The RS model is a 5-dimensional theory with the fifth dimension compactified on 
an /Z2 orbifold with a radius Rc- The Planck brane is located at the orbifold fixed 
point (f) = 7T while the SM fields are localized at the TeV brane which is at = 0. This 
geometry gives the following metric in 5-dimensions: 

where < < vr. To explain the hierarchy between the Planck scale and the electroweak 
(EW) scale we need KRc only of the order (9(10). Introducing an extra scalar field in 
the bulk [11[5] showed that JCRc can be made stable against the quantum fluctuations. 

The variations of the above setup have also been considered in the literature where 
the SM fields, except for the Higgs field, have been allowed to propagate in the bulk P-[8] . 
This framework provides an interesting new approach to the flavor problem, as now also 
the hierarchical structures observed in the masses and the mixing of the SM fermions 
could be explained in terms of geometrical effects [7] , [814TT] . We will consider the original 
proposal of RS for our analysis. 

The effect of extra dimensions on the SM fields is felt through the KK gravitons. 
These KK gravitons, h^^J , couple to the SM energy momentum tensor and the interaction 
Lagrangian is 

Mp, ^ ' t^''^ ' Mp, ^ \ J J \ J 
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T'^'^ is the symmetric energy-momentum tensor for the SM particles on the 3-brane, and 
Mpi is the reduced Planck scale. The masses of the h^^J are given by 

M„ = XnIC e--'^«= , (3) 

where the x„ are the zeros of the Bessel function Ji{x). The first term in the interaction 
Lagrangian gives the coupling of the zero-mode and it is Planck scale suppressed. The 
coupling of the massive KK states is enhanced due to the exponential factor e'"'^^" and 
gives interactions of EW strength. Consequently, except for the overall warp factor in 
the RS case, the Feynman rules in the RS model are the same as those for the ADD 
case [T2l[T3l[27] . The basic parameters of the RS model are 



mo = ICe 



-TtKRc 



Co = IC/Mpi , (4) 

where mo is a scale of the dimension of mass and Cq (0.01 < Cq < 0.1) is an effective 
coupling. For our analysis we choose to work with the RS parameters cq and Mi the 
first excited mode of the graviton rather than mo. 

Summing over all the KK states we obtain effective graviton propagator : 

= ^ Q^^^M^TO^ ^ ^ ' 

where M„ are the masses of the individual resonances (see Eq. [3]) and the r„ are the 
corresponding widths. 

There are two ways to probe such effects at colliders, either through graviton emission 
or by virtual graviton exchange. In this paper we will consider only virtual spin-2 KK 
states. Production of boson pairs is one of the important process at the LHC both in the 
context of SM and new physics studies. Studies in other channels have been reported in 
[H] in extra dimension models. In this paper we will consider production of W pair at the 
LHC. Owing to its importance many studies have been carried out for its production in 
the SM; a study in the context of anomalous triple gauge boson vertices was carried out in 
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p3|[T6]. Leading order (LO) studies in the SM can be found in [T7]. As is well known the 
LO results are highly sensitive to the arbitrary renormalization and factorization scales. 
At this order the factorization scale fip enters solely through the parton distribution 
functions as the parton level cross-section, at this order, does not depend on fip. As we 
include higher order terms of the perturbation series the dependence will reduce and an 
all order result will be completely independent of these arbitrary scales. In addition the 
NLO results are usually significantly enhanced as compared to the LO results. It is thus 
important to carry out NLO calculation to reduce these scale dependencies. Because of 
its importance, its production has been studied to next-to-leading-order (NLO) accuracy 
in the SM [T81 - I2T] . These results were subsequently updated in [22l[23]. These studies 
provide the precise estimate of higher order effects through K factor as well as the 
sensitivity of the predictions to factorization scale. Its production has also been studied 
via gluon fusion through a quark box loop or triangle quark loop with 7 or Z boson 
exchange |23] and at one and two loop level in high energy limit in SM |25] . 

Two W bosons can couple to Kaluza Klein {KK) gravitons, so it is possible to 
produce them through virtual graviton exchange at LO [26j. The significance of NLO 
computations in the extra dimension models for Drell-Yan [27], diphoton [28], ZZ [29] . 
graviton+photon [30], graviton -(-jet [31] production has already been demonstrated. Al- 
though NLO results are available in SM, they do not exist in literature in the context of 
RS model for W boson pair production, which is the material of the present paper. 

Before we present the results let us present in brief the pieces of NLO calculation. 
The details can be found in [32] where we have given the matrix elements etc. for the 
process in the context of ADD model. The signal comprises of contributions 

\Msm\^ + I-MgI' + iMsMM*a + c.c.) (6) 

where the first term is pure SM, the second is purely gravity mediated and the third term 
is the interference of SM and gravity mediated processes. At leading order in strong 
coupling AisM has three contributions; a t-channel or u-channel process and s-channel 
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processes via 7 and Z boson. 

qq'^W^W, qq'AW^W-, qq "4 W+W- (7) 

As the KK gravitons couple with same strength to quarks and gluons both quark and 
gluon initiated Feynman diagrams with s-channel graviton propagator contribute to A^g- 

Next at order as we have to include both one loop corrections to the above processes 
and also real emission contributions in which in addition to W~^W~ a parton is emitted in 
the final state. The soft and collinear configurations in the loop integrals give divergences 
which we have regulated using dimensional regularization [n = 4+e) thus the singularities 
appear as simple and double poles in e. As the process under consideration is UV finite, 
these poles are only soft and collinear. In the real emission case we have qq, qg and 
gg initiated processes. As we have gg — > W'^W" at leading order through graviton 
exchange, we note that all the 4-kinds of splitting functions Pgg, Pgg, Pgg, Pgg are involved 
in the calculation. In addition to the above soft and collinear singularities, the other 
set of these divergences appear from phase space integration of the real emission matrix 
elements. The sum of virtual and real contributions is completely finite ie. free of poles 
in e after mass factorization is carried out. We have used MS scheme throughout, both 
for the renormalization and factorization. 

We have employed the method of two cutoff phase space slicing to handle the real 
emission processes. In this method two small dimensionless slicing parameters 5s and 5c 
are introduced to divide the real emission phase space into soft and collinear regions (for 
a review of the method please see 03]). The cross section can be written as, then, 

da = da^° + rfa"''* + 4) + da"^""'^ """^ 4) (8) 

Here the third term gives the contribution coming from the soft and collinear regions 
which is rendered finite after adding the counter term (CT) for mass factorization. The 
last term denotes the contribution of hard non collinear configurations and is finite. We 
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define 



(9) 



hard non col 



(10) 



Note that, individually rfa^"''"'^^ and c/cr^"'""^^ depend on 5^ and Sc but the sum should 
be independent of the parameters which were introduced to slice the phase space. We 
have incorporated all the above details in our monte carlo code which is implemented on 
FORTRAN 77 and easy to tailor for various cuts on the final state bosons. 

We now make some general remarks about the computation. We have used Feynman 
gauge in QCD sector and unitary gauge in electroweak sector. The choice of unitary 
gauge simphfies the calculation as both the electroweak Goldstone bosons and ghosts 
disappear. Further we note that the term proportional to 1/.^ in gluon-gluon-graviton 
vertex does not contribute. Also the results do no depend on the arbitrary vector n^^ 
which appears in gluon polarization sum: 



Further our SM matrix elements agree with those given in p^|[T9] . To check the numerical 
implementation of the phase space slicing method we have checked the stability of the 
sum of 2 — body and 3 — body contributions against variation of slicing parameters 6s and 
6c and we found the sum to be stable over a wide range of these parameters. In what 
follows we will use 6s = 10~^ and 6c = 10~^. 

We now present the kinematical distributions for the W~^W~ production at the LHC 
The LHC with a center of mass energy of 14 TeV will be our default choice. However 
we will also present some results for a center of mass energy of 10 TeV for the LHC. For 
numerical evaluation, the following SM parameters [31] will be used 



mw = 80.398 GeV, mz = 91.1876 GeV, Tz = 2.4952 GeV, sin^ Ow = 0.231 (12) 



where 6\y is the weak mixing angle. For the electromagnetic coupling constant a we 
use = 128.89. CTEQ6 density sets are used for parton distribution functions. 



k ■ n 



(11) 
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2-loop running for the strong coupling constant is used . The number of active hght- 
quark flavors is taken to be 5 and the value of Aqqd is chosen as prescribed by the 
CTEQ6 density sets. At leading order we use CTEQ6L1 density set ( which uses the 
LO running ) with the corresponding Aqcd = 165 MeV. At NLO we use CTEQ6M 
density set ( which uses 2-loop running ) with the Aqcd = 226 MeV; this value of 
Aqcd enters into the evaluation of the 2-loop strong coupling. The default choice for 
the renormalization and factorization scale is the identification to the invariant mass of 
the W boson pair ie., fip = fiR = Q. Furthermore the W bosons will be constrained to 
satisfy \yw\ < 2.5, where yw is the rapidity of a final state W boson . 

We present invariant mass (Q) and rapidity (Y) distribution of the W boson pairs. 
These kinematical variables are defined as 

Q' = iPw^+Pw-r, Y = \\n^, (13) 

where Pi and P2 are the momenta of colliding hadrons, and q = pw+ + Pw- denotes the 
sum of the VT-boson 4-momenta. In obtaining these distributions all order as contribu- 
tions have been taken into account. 

In Fig. [T] we have plotted the invariant mass distribution both for the SM and the 
signal for LHC at lATeV. The two curves with peaks correspond to the signal and 
the remaining two curves give SM predictions. Here we have chosen cq = 0.01 and 
Ml = 1500GeV. To highlight the importance of QCD corrections we have also displayed 
the LO results of SM and the signal, and we observe that at Q = 1500 GeV the K 
factors (defined as i^' = da^^^ / da^^) has a value 1.9. Thus NLO QCD corrections give 
a substantial enhancement over the LO predictions. 

Next we present in Fig. [2] the effects of varying the parameter cq on the invariant 
mass distribution. All the curves shown correspond to NLO results, and we have also 
plotted the SM background for comparison. 

In Fig. Owe have plotted the rapidity distribution da/dY at NLO both for SM and 
the signal for cq = 0.01. We have plotted this distribution in the interval —2.0 < Y < 2.0 
and have carried out an integration over the invariant mass interval 1450 < Q < 1550 to 
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increase the signal over the SM background. As expected the distribution is symmetric 
about Y = 0. 

As was noted above the NLO QCD corrections reduce the sensitivity of the cross 
sections to the factorization scale fip; this we now show in the Fig. HI We have plotted 
SM and the signal both at LO and NLO, and have varied the factorization scale fip in the 
range Q/2 < fip < 2Q. The central curve in a given band (shown by the dotted curves) 
correspond to fip = Q. In all these results the renormalization scale is fixed at = Q. 
We notice that the factorization scale uncertainty at LO is 21.8 % at Q = 1500 GeV 
as compared to 6.7 % at NLO. Thus we see that NLO computation achieves significant 
reduction in uncertainty and makes predictions much more precise. 

At the end we present in Fig. [5l da/dQ for LHC with a centre of mass energy of 
10 TeV at NLO both for SM and signal. For comparison we have also plotted the 14 TeV 
results in the same figure. 

To summarize, in this paper we have carried out a full NLO QCD calculation for the 
production of two W bosons at the LHC at 14 TeV and 10 TeV in the extra dimension 
model of Randall and Sundrum. Here we take all order contributions, both in the 
SM and in the gravity mediated processes and their interferences, into account. We 
have presented invariant mass and rapidity distributions both at LO and NLO. We use 
CTEQ 6L1 and CTEQ 6M parton density sets for LO and NLO observables, respectively. 
Significant enhancements over the LO predictions are observed. The K factor are large 
and at Q = 1500 GeV (we have taken this as the first RS resonance) K = 1.9. This 
justifies the entire exercise of carrying out a NLO computation. The effect of variation 
of parameter cq in invariant mass distribution is also presented. We have shown that a 
significant reduction in LO theoretical uncertainty, arising from the factorization scale, 
is achieved by our NLO computation. It is observed that an uncertainty of 21.8 % at 
LO as /ii? is varied between Q/2 and 2Q is reduced to 6.7 %. Thus our NLO results are 
more precise than the LO results and suitable for further studies for constraining the 
parameters of the RS model. Invariant mass distribution is also presented for LHC at a 
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center of mass energy of lOTeV at the NLO level. 
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Figure 1: Invariant mass distribution for SM and signal both at LO and NLO. Dash- 
dot curves represent LO results and solid curves give NLO results. We have chosen 
Ml = 1500 GeV and the parameter cq = 0.01. 
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Covariation, LHC (V S =14 TeV) 
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Figure 2: Effect of variation of cq on invariant mass distribution. All the curves corre- 
spond to NLO results with Mi fixed at 1500 GeV. The solid curve corresponds to SM 
and the dash-dot curves to the signal. The signal is plotted for cq = 0.01, 0.04, 0.08 and 
the dash size increases with increasing cq 
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LHC (V S =14 TeV) 
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Figure 3: Rapidity distribution for SM and signal at NLO. Dash curve represents the 
signal and solid curve gives SM result. We have chosen Mi = 1500 GeV and the pa- 
rameter Co = 0.01. To enhance the signal we have integrated over Q in the range 
1450 <Q< 1550. 
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Hp variation, LHC (V S =14 TeV) 




Figure 4: Factorization scale variation in the invariant mass distribution. The curves 
correspond to Cq = 0.01 and Mi = ISOOGeV" at the LHC at = 14 TeV. The is 
varied between Q/2 and 2Q. The dash curves correspond to /if = Q 
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LHC (V S =14 TeV, 10 TeV) 
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Figure 5: Invariant mass distribution for SM and signal at y/S = lOTeV and lATeV. 
All the curves correspond to NLO results. We have chosen Mi = 1500 GeV and the 
parameter cq = 0.01. 
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